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Electromagnetic methods recently proposed for detecting gravitational waves modify the Michel- 
son phase shift analysis (historically employed for special relativity). We suggest that a frequency 
modulation analysis is more suited to general relativity. An incident photon in the presence of a very 
long wavelength gravitational wave will have a finite probability of being returned as a final photon 
with a frequency shift whose magnitude is equal to the gravitational wave frequency. The effect 
is due to the non-linear coupling between electromagnetic and gravitational waves. The frequency 
modulation is derived directly from the Maxwell-Einstein equations. 
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I. INTRODUCTION 

There has been considerable recent interest in the na- 
ture of gravitational and electromagnetic wave interac- 
tions, especially regarding optical techniques for detect- 
ing gravitational waves [lflllSllSIEHH. Among the 
possible electromagnetic techniques for detecting gravi- 
tational waves is a modification of the methods which 
were historically important for the application of Michel- 
son interferometers |9j to special relativity. A new tech- 
nique can be employed as central for the detection of 
gravitational waves. It resides in the non- linear coupling 
between an incident traveling gravitational wave and al- 
most standing Fabry-Perot cavity electromagnetic waves. 
Theoretically, such a non-linear wave coupling leads to a 
dynamic modulation and frequency shift of a finally de- 
tected electromagnetic signal. Such a dynamic (for gen- 
eral relativity) frequency shift measurement would super- 
sede the earlier notions of Michelson who concentrated 
on static (for special relativity) phase shifts as measured 
from an interference pattern. Recall the Michelson inter- 
ferometer as schematically pictured in Fig^ 

Our purpose is to provide a detailed derivation of the 
gravitational modulation of the electromagnetic cavity 
waves from the coupled Maxwell-Einstein theory. In 
Sec[n] we derive the tree level Feynman diagram corre- 
sponding to a photon transition 7, — > 7^ in the presence 
of a gravitational wave. If, respectively, u> g (u>i) denotes 
the frequency of the gravitational wave (initial photon), 
then the final photon frequency Wf in the modulation 
sideband obeys 



UJf = UJi ± LUg 



(1) 



In Sec lIIII the electromagnetic modes of the cavity in the 
presence of a gravitational wave are explored. The fre- 
quency modulation follows from the nature of the cou- 
pling between the gravitational wave and the Maxwell 
electromagnetic pressure tensor. In Sec II VI the "Fermi- 
Golden- Rule" rates for the transitions 7, — > 7/ in Eq.QJ 
are computed in terms of the power spectrum of grav- 
itational strain fluctuations, and both electromagnetic 
and gravitational polarizations will be discussed. In the 



concluding Sec0 the notion of gravitational wave in- 
duced modulation side bands will be discussed in virtue 
of the vacuum polarization response induced by the grav- 
itational strain. 



II. ELECTROMAGNETIC AND 
GRAVITATIONAL WAVE INTERACTIONS 

The mathematical form of the electromagnetic- 
gravitational wave interaction arising from the gravita- 
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FIG. 1: Shown schematically (and not to scale) is a modern 
day Michelson interferometer built for gravitational wave de- 
tection. One seeks to measure the spatial strain of an incident 
gravitational wave whose wavelength X g — (2irc/oj g ) is large 
on the scale of the apparatus size. An incident photon 7^ can 
take one of two distinct paths into (and out of) the vertical 
or horizontal Fabry-Perot cavities. The final detected photon 
7/ can (with finite probability) have its frequency shifted into 
a side band by the amount determined by the gravitational 
wave frequency; i.e. uif — UJi = ±ui g . 
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tional metric 

— c 2 dr 2 = g^dx^dx" , 
and the cavity electromagnetic field 

Fp,v = d^An — d u A^ n 
can be described by the action 



(2) 
(3) 

(4) 



wherein dtt = y/—g d^x. For a small change in the grav- 
itational metric 8g^ v representing a weak incident gravi- 
tational wave, the coupling into the electromagnetic field 
may be written in terms of the Maxwell stress tensor T MI/ ; 



(5) 



Since the Maxwell stress tensor is quadratic in the elec- 
tromagnetic fields, Eqs.|J5J determine the wave coupling 
strengths between the two photons (initial and final) and 
a gravitational wave in accordance with the Feynman di- 
agram of Fig|21 

During the time period in which the photons are part 
of an almost standing wave within a Fabry-Perot cav- 
ity pictured in FigQ] a scattering event may take place 
where the initial photon either absorbs (from the grav- 
itational wave) or emits (into the gravitational wave) a 
single quantum Tioj g of energy. Such gravitational wave 
emission or absorption processes will yield photon side- 
bands of frequency, Wf — ioi ± uj g . Since the Feynman 
diagram in Figgis a "tree" (with no internal "loops"), 
it follows that a classical non-linear wave approach to 
the problem yields answers equivalent to the more fun- 
damental quantum theory. 



So 



FIG. 2: Shown is a process 7, 7/ in the presence of an 
external gravitational metric disturbance <5g M „. The Feynman 
diagram describes the scattering of light in a gravitational 
field. For a gravitational wave disturbance at frequency u) g , 
the frequency shift in the photon obeys Uf — Wj = ±w g . Such 
a frequency shift should appear experimentally as modulation 
sidebands about the initial photon frequency Wj. 



III. FREQUENCY MODULATION 

For the purpose of describing the gravitational wave, 
we employ the metric expansion in which g^ differs only 
slightly from flat space-time; i.e. 



9nv V^i-v "h h[n/ ~\- 



(6) 



where flat space time is described by rj^ and is de- 
scribed by a spatial transverse traceless strain u; i.e. 
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(7) 



From Eqs.Q and JSJ it follows (to lowest order in the 
gravitational strain u) that the Maxwell field Lagrangian 
is 



LMaxwell — ^— I (|E| 



(8) 



where the interaction between the electromagnetic field 
and the gravitational strain is given by 



L mt = J(u : P)d 3 r. 



(9) 



In Eq.©, P y = T" y is the electromagnetic spatial pres- 
sure tensor, normally written as 

P= i-{l(|E| 2 + |B| 2 ) -2(EE + BB)}. (10) 

Since we work in a traceless gauge tr(u) = 0, the La- 
grangian coupling of the Fabry-Perot cavity fields (E and 
B) to the gravitational strain u is simply written as 



Li 



J (E • u • E + B • u ■ B)d 3 r. (11) 



We note that an "effective" electromagnetic La- 
grangian L[E, B] may be used to define the Maxwell dis- 
placement field D and the magnetic intensity H via the 
functional derivatives 



SL 

D = An— and H 
<5E 



-47T- 



SL 
SB 



(12) 



Thus, Eqs.JHJ, (fTTJl and (ft2"|l imply a tensor dielectric re- 
sponse D = e ■ E and a tensor magnetic permeability 
B = /i ■ H determined by the gravitational wave; Explic- 
itly, the response functions are given to linear order in u 
by 



e jk (r,t) = 6jk - 2u jk (r,t), 
(^ 1 ) ]k {T,t) = 5jk+2ujk(r,t). 



(13) 



A gravitational wave thereby acts (via e and fx) as a 
weak "moving grating" in the vacuum modulating the 
frequency of traveling electromagnetic waves. 
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In the limit of long wavelength gravitational waves, the 
strain u is uniform in space over the length scale of the 
interferometer. One thereby can apply the gravitational 
wave quadrupole approximation and consider that u de- 
pends only on time. Eq.© then reads 



L int = u(t) : J P(r,i)d 3 r. 



(14) 



To lowest order perturbation theory in u, the interac- 
tion Lagrangian in Ea. (|14l) may be replaced by the 
Hamiltonianlld 



H^t = -u(i) : J P(r,i)d 3 r. 



(15) 



It is theoretically fortunate (and experimentally unfortu- 
nate) that the gravitational wave induced photon transi- 
tions are so weak that the lowest order in u perturbation 
theory is virtually exact. 



IV. PHOTON TRANSITION RATES 

The amplitude for a transition to take place in a long 
time period r is given in lowest order perturbation theory 
as 



A(I^F;t) 



(F\H mt (t)\I) alt 



u(t) : Q F ie iUFlt dt, (16) 



where the Bohr transition frequency ujpj = (Ep — Ex)/h 
and 



Qfi = J (F\ P(r) | J) d 3 r. 



(17) 



The transition rate per unit time for a radiation tran- 
sition, 



r 



I—>F 



lim 



\A(I~>F;t) 



follows from Eas. (|16[) and (|18fl to have the form 
2tt 



I—>F 



h 2 



(Qjk) FI (Qlm) pjSjklm {UFI 



(18) 



(19) 



Apart from the matrix elements squared of the pho- 
ton field operators / in Eas. (|17(l . the effective (fourth 
rank tensor) "density of states" in the Fermi golden rule 
Eq. (|19JI is determined by the power spectrum of gravita- 
tional wave strains via 



Sjkim{u) = lirrir^oo (u jk (uj)*ui m (w)) (20) 



where 



u(uj) 



e iuJt u(t)dt. 



(21) 



The power spectrum may equally well be written in terms 
of the two time correlation functions 



«im(*a)«i*(*i) = / S jk i m (oj)e- iu ^-^(kj. (22) 



Let £ be the polarization tensor for a gravitational wave 
described by the strain u. Furthermore let 



(23) 



be the power spectrum of gravitational strain fluctuations 
with polarization £. Such a power spectrum is easily 
related to the gravitational wave energy per unit time 
per unit area dV^(oj) incident on the interferometer in a 
bandwidth du>; It is (for u> > 0) 



dV^u 



\4nG 



S^(u>)du>, 



(24) 



where G is the gravitational coupling strength. 

For a gravitational wave emission source with a fre- 
quency Lo g and a Lorentzian width 7, the resulting power 
spectrum is plotted in FigEI The incident gravitational 
wave acts as an effective density of final states for the 
induced photon transition rate 7$ — > 7/. If the gravita- 
tional wave intensity distribution (dV^(uj) / dm) is singly 
peaked, then the "density of final states" will have two 
symmetrically located peaks. These will give rise to the 
modulation sidebands in the spectrum of detected pho- 
tons emerging from the interferometer. 

The Fabry-Perot detector efficiency may be computed 
from Ea. H19(l by averaging over initial radiation states 



Modulation Side Bands 
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FIG. 3: Shown is the power spectrum S^(u>) of gravitational 
strain oscillations corresponding to a frequency ui g and a 



Lorentzian emission width 7. 
oscillation is u 2 and (wg/7) = 
trative purposes. 



The integrated strength of the 
= 3.0 has been chosen for illus- 
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and summing over final radiation states; i.e. 



V. CONCLUSION 



dT(u) 
dT(u) 

did 



I F 



(25) 



where the radiation (detector) correlation functions are 
defined as 



Gjklm (t) 
Qjklm{tl — t\) 



V jklm (Lo)e- luJt dw, 



% (Qjk(ti)Qi m {t 2 )) (26) 



Since Q(t) = J P(r,t)d 3 r, the detection efficiency is, in 
virtue of Eas. l|25|l and i|26[) . proportional to the power 
spectrum of the integrated radiation pressure in the 
Fabry-Perot cavities. If the bandwidth of the input laser 
is large on the scale of the gravitational wave frequency, 
then the detection efficiency is roughly frequency inde- 
pendent and oc N^/Hz where iV 7 is the number of pho- 
tons stored in the Fabry-Perot cavities. The signal to 
noise ratio for detecting the modulation side bands is 
given by 



signal y AirS^ (a 



noise 



(27) 



where ~ 10~ 22 / V Hz is the strain noise level with 
present technologies for kilohertz gravitational wave sig- 
nals. 



We have shown that the electromagnetic waves in the 
Fabry-Perot cavities of a Michelson interferometer may 
be described by the Lagrangian density 



f 



£= — (E-e-E-B-/* _1 -B) 



(28) 



wherein the vacuum dielectric response and magnetic 
permeability are related to the gravitational wave trans- 
verse traceless tensor via Ea. lfFl)) : i.e. 



(/i^jfcfr,*) = S jk + hjk(r,t). 



(29) 



The coordinates r being used are such that the ends of the 
cavities have fixed positions in a gravitational "floating" 
situation. Since the e and fi tensors formally change light 
velocity in the arms of the interferometer, the distances 
between the ends of the arms (as measured by travel 
times of light signals) will be changing when a gravita- 
tional wave is present. An oscillating e and /i will then 
produce frequency modulation side bands. Such modula- 
tions constitute a very well known effect [TlT| in quantum 
optics. 
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